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C"| ■ Abstract. We explore squeezed coherent states of a 3-dimensional generalized 

O-i isotonic oscillator whose radial part is the newly introduced generalized isotonic 

oscillator whose bound state solutions have been shown to admit the recently discovered 
Xi-Lagucrrc polynomials. We construct a complete set of squeezed coherent states 
^ | of this oscillator by exploring the squeezed coherent states of the radial part and 

O"* combining the latter with the squeezed coherent states of the angular part. We also 

prove that the three mode squeezed coherent states resolve the identity operator. We 
evaluate Mandel's Q-parameter of the obtained states and demonstrate that these 
states exhibit sub-Possionian and super-Possionian photon statistics. Further, we 
qq | illustrate the squeezing properties of these states, both in the radial and angular parts, 

by choosing appropriate observables in the respective parts. We also evaluate Wigner 
function of these three mode squeezed coherent states and demonstrate squeezing 
property explicitly. 



PACS numbers: 03.65.-w, 03.65. Ge, 03.65.Fd 
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1. Introduction 

1.1. New exactly solvable model 

Very recently the exact quantum solvability of the extended radial oscillator/ generalized 
isotonic oscillator potential, 

, 22 1(1 + 1) 8cu 16u(2l + 1) , 1X 

V(r) = cu r + — H — , ( 1) 

{} r 2 (2wr 2 + 2/ + l) [2ur 2 + 21 + l) 2 ' K) 

where u and / are parameters, has been demonstrated in four different perspectives 
[DEI El II]- To begin with, it was shown that the bound state eigenfunctions of the 
Schrodinger equation associated with the potential ([T|) can be expressed in terms of 
the newly found exceptional orthogonal v th degree polynomials, L k (x), v — 0, 1, 2, 3, ... 
where A; is a positive real parameter, namely Xi-Laguerre polynomials [I]. These new 
polynomials, L k (x), are related to the classical Laguerre polynomials L k (x) by the 
following relation [5] 

l\(x) = -(x + k + l)Ll_ x {x) + L k v _ 2 (x). (2) 

The bound state solutions of the Schrodinger equation associated with the potential 
(CQ) are found to be [1] 



E n j = 2u ( 2n + l + - ) , n = 0,1,2,3..., I = 0,1,2,3, ... (4) 



with 

E„ , = 2ui ( 2n 4- 1 4- ^ 

1/2 

and A^; = I — ^ ) is the normalization constant. We note here 

{n + l + |)r(n + I + ± ; 

that the Xi-Laguerre polynomials, L k (x), are solutions of the following second-order 
linear ordinary differential equation with rational coefficients, that is 

^„ _ (x-k)(x + k + l) , + 1 / x - k + v _^\ Q ' = A. (5) 

x(x + k) x \x + k ) dx 1 

where k > is a real parameter and v = 1,2,3,... . The first few XpLaguerre 
polynomials are [5] 

L k {x) = -(x + k + 1), (6) 
L k (x) = x 2 -k{k + 2), (7) 

L k (x) =-\x* + + _ * (3 + 4k + ( 8 ) 

One may observe that the above polynomial sequence starts with a linear polynomial 
in x instead of a constant which usually the other classical orthogonal polynomials 
do. However, these new polynomials form complete set with respect to some possible 
measure [5j. For more details on the properties of these exceptional orthogonal 
polynomials, one may refer Ref. [5]. 
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In a different context, while constructing exact analytic solutions of the d- 
dimensional Schrodinger equation associated with the generalized quantum isotonic 
nonlinear oscillator potential ([T|), that is 

where A is the d-dimensional (d > 2) Laplacian operator and u, a, g are parameters and 
V(r) is the central potential given by 

V{r) = ^ + ^ + 2g f- a l ) B 2 >0, (10) 

Hall et al have shown that only for the values g = 2 and u 2 a 4 = B 2 +(l + ^) the 
potential ({TO]) can be exactly solvable |2J. It is shown that the potential ( !T0|) can be 
regarded as a supersymmetric partner of the Goldmann - Krivchenkov potential whose 
exact solutions are known [2]. The potential f JTQj) with the parametric restrictions given 
above exactly matches with the one constructed by Quesne [JJ. In a subsequent study [3J, 
Saad et al have solved the system QUI) numerically, for arbitrary values of the system 
parameters, through asymptotic iteration method. They have shown that for certain 
restrictions (B — and d — 1) the potential (flOl can be reduced to the one discussed in 
Ref. [6J (see equation ffTTj) below and the paragraph that follows). The quasi polynomial 
solutions and the energy eigenvalues of ( TTOl are also reported in Ref. [3j. 

In a very recent paper, Agboola and Zhang have considered the potential ([1]) and 
transformed the associated Schrodinger equation into the generalized spheroidal wave 
equation [7]. Using Bethe ansatz method, they have derived bound state solutions of 
(JTJ) for certain parametric values jl]. 

We mention here that the inverse square type potentials play an important role in 
physics [8|l9t lTmfTTp2] . Solving Schrodinger equation with non-polynomial potentials [13] 
is also of interest for physicists at both the classical [H] and quantum level [T5] . 

We also recall here that during the past few years the one dimensional version of 
the potential (JTJ) , that is 

7(I)=W¥ + 2 ^' a2> °' (U) 

is studied in different perspectives, see for example Refs. [3j[6l[T^[T^l[T8l[T9l|20l|2T]. The 
generalized isotonic oscillator potential (II ip is generalized in such a way that it lies 
between the harmonic oscillator and isotonic oscillator potentials [6J. Importantly, this 
generalization removes the singularity nature of the isotonic oscillator in real space, 
as the generalized isotonic oscillator has poles only at imaginary points x = ±ia. 
The energy values of this generalized isotonic oscillator potential (g = 2ua 2 (l + 
2ua 2 ), u = 1, a 2 = \) are equidistant [B]. In our earlier studies, we have constructed 
different types of nonlinear coherent states, nonlinear squeezed states and various 
nonclassical states of the oscillator ( II ip and analyzed its position dependent mass 
Schrodinger equation by fixing the parameters g = 2ua 2 (l + 2ua 2 ), u — 1 and 
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a 2 = | [201 1211 E21 EH]- While constructing nonlinear coherent and squeezed states of 
the system (ITT]), the deformed ladder operators which we had constructed lead to two 
non-unitary displacement operators and two squeezing operators [2TJ [22l |23] . While 
one of the displacement operators produced nonlinear coherent states the other non- 
unitary displacement operator failed to produce any new type of nonlinear coherent 
states (dual pair) [23] . In the case of squeezed states also, while one of the squeezing 
operators gave the nonlinear squeezed states the other one failed to produce the dual 
pair of nonlinear squeezed states which lead us to the conclusion that the dual pair of 
nonlinear coherent and squeezed states are absent in this system (for more details one 
may refer the Refs. (22lj23]). 

1.2. Present work 

As a continuation of our earlier studies, in this paper, we focus our attention on the three 
dimensional extension of the potential ( ITT]) . With the restriction d = 3, g = 2, B = 
and coa 2 — I + |, the Schrodinger equation ([9]) turns out to be an exactly solvable three 
dimensional isotonic oscillator equation, that is 

- + u 2 r 2 + ^ — - — '— = EV. 12 

V 2wr 2 + 21 + 1 (2ur 2 + 21 + I) 2 ) K J 

The isotonic oscillator with the generalized inverse type of potential models (i) the 
dynamics of harmonic oscillator in the presence of dipoles where it describes the 
scattering off electrons by polar molecules [25J and (ii) the dynamics of dipoles in Cosmic 



strings [TT] . 

Now considering the three dimensional Laplacian in polar coordinates with the 
wavefunction \l/(r, 6,4>) = u{r) Y^ m (9,<p), where Yi tTn (6,4>) are the spherical harmonics, 
we have the Schrodinger equation 

d 2 u _2du f 22 1(1 + 1) 8u _ 16w(2/ + l) \ 

~ d^ 2 ~ rdr~ + \ T + r 2 + 2ur 2 + 21 + 1 ~~ {2ur 2 + 2l + l) 2 ) U ~ ( } 
where / is the angular momentum. 

Since the bound state solutions of (|T3|) are already known (vide equations (j3J) and 
(SJ)) we can unambiguously write the eigenf unctions and energy values of the three 
dimensional problem as [1,2J 



<Mr,0,# = N nil — 2 ^ 2f + 1 4+i*W)e-* w Xm(M), (14) 



E n> i = 2cu \ 2n + I + , n = 0, 1, 2, 3, 1 = 0, 1, 2, 3, m = -I to I. (15) 

From these solutions we construct a complete set of squeezed coherent states for the 
three dimensional system ( fl2l) . Since the potential under consideration is a spherically 
symmetric one, we separate radial part from the angular part and investigate each one of 
them separately. To construct squeezed coherent states of the radial part, vide equation 
([I]), we should know the ladder operators of it. The ladder operators can be explored 
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through supersymmetric technique. However, the resultant operators which come out 
through this formalism do not perform perfect annihilation and creation actions on 
the number states. To explore the perfect ladder operators we alternatively examine 
shape invariance property of the potential ([!])■ Once the presence of this property 
is confirmed, we then move on to construct the necessary annihilation and creation 
operators by adopting the procedure described in Ref. [26]. Using these operators 
we unambiguously construct squeezed coherent states [271 [28] of the radial part. To 
explore squeezed coherent states of the angular part we use Schwinger representation and 
define the necessary creation and annihilation operators. In particular, by considering 
disentangled form of two mode squeezing and displacement operators we derive squeezed 
coherent states [29]. We obtain a complete set of squeezed coherent states of the 3- 
dimensional oscillator by defining these states are the tensor product between squeezed 
coherent states of the radial variable and squeezed states of the angular variables [30J. 
To make our result more rigorous we also prove that these three mode squeezed 
coherent states satisfy the completeness condition. We also evaluate Mandel's Q- 
parameter of the three mode squeezed coherent states and investigate certain photon 
statistical properties associated with these states. Our result confirms that these states 
exhibit sub-Poissonian (non-classical) and super-Poissonian photon statistics for certain 
parameters. Further, we demonstrate squeezing properties of the states associated with 
the radial part by considering the generalized position and momentum coordinates 
with the obtained ladder operators and the angular part by considering the angular 
momentum quantities L x and L y respectively. Finally, we evaluate the Wigner function 
of the constructed squeezed coherent states by defining a three-mode Wigner function in 
terms of coherent states [HUGS]. Our result reveals that these squeezed coherent states 
exhibit the squeezing property. We intend to carry out all these studies due to the 
fact that among the nonclassical states squeezed states have attracted much attention 
because certain observables in this basis show the fluctuations less than that of the 
vacuum. The squeezed coherent states have found applications not only in the quantum 
optics but also in quantum cryptography [33J, quantum teleportation [31] and quantum 
communication |35j, to name a few. 

This paper is organized as follows. In section 2, we analyze the shape invariance 
property of the radial part of the 3-dimensional generalized isotonic oscillator and 
construct suitable ladder operators. In section 3, we construct squeezed coherent states 
of the three dimensional oscillator. We use these ladder operators to construct squeezed 
coherent states of the radial part and Schwinger's operators to construct angular part. 
We express the complete form of the squeezed coherent states of the 3-dimensional 
oscillator by suitably combining the radial part with the angular part. We also prove 
that the three mode squeezed coherent states satisfy the completeness condition. We 
investigate certain non-classical properties associated with these states in section 4. To 
begin with, we study Mandel's Q-parameter and confirm the non-classical and super- 
Poissonian nature of these states. We also illustrate the squeezing properties of the 
constructed 3-dimensional squeezed coherent states both in the radial part and in the 
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angular part in this section. In section 5, we evaluate Wigner function for the 3- 
dimensional generalized isotonic oscillator and demonstrate the squeezing property of 
the constructed states. Finally, in section 6, we present the conclusion and the outcome 
of our study. Certain background derivations which are needed to obtain the necessary 
ladder operators are given in the Appendix A. The details of evaluating certain multiple 
sums in the expressions (n r ) and (n 2 .) are given in Appendix B. In Appendix C we 
discuss the method of evaluating the multiple sums appearing in the Wigner function. 



2. Shape invariance property and ladder operators 

2.1. Supersymmetric formalism 

To construct annihilation and creation operators, we start with time-independent 
Schrodinger equation, 

H^r) = (^-— + V{r) S j<$>{r) = E®{r), (16) 

with V(r) is given in (CQ). Let us rewrite the second order differential operator H given 
in fll6p as a product of two first order differential operators, namely A + and A~, such 
that #W = H - u(2l + 3) = A+A~ and 

A + = -— + W(r) and A~ = ^- + W(r). (17) 
dr dr 

In the above W(r) is the superpotential which is found to be [TJ 

■art \ 1 + 1 i 4wr 4wr fia\ 
W(r)=ur . (18) 

v ; r 2ur 2 + 21 + 1 2ur 2 + 21 + 3 ( J 

The operators A + and A~ factorize the partner Hamiltonian H^ 1 as 

H {2) = A~A + = H-u{2l + \). 

The Schrodinger equation associated with reads 

ff«*g!(r) = £f«$gi(r) > (19) 

that is 

where the potential V\{r) is given by 

— — , 22 1(1 + 1) 8u 16u(2l + l) , nt oN . . 

Equation ( 120]) shares the same solution as that of the Schrodinger equation associated 
with the potential (OQ) admits, that is 

*S s = N -' ^/+2l + l) L ^ ] ^ ^ (22) 
with the energy eigenvalues 

= E n;l -u(2l + 3)= Anu. (23) 



rvu/p ni ,v;-^'* nil \f), (20) 



The Schrodinger equation associated with is given by 

^ (2) < 2 lW = ^ ) $S(r). 

Rewriting (124"]) . we get 



dr 2 



where the partner potential V^r) turns out to be 

W N 2 2, (/ + !)(/ + 2) , 8w 
(r) = + - 



16w(2/ + 3) 



(24) 
(25) 

u(2l + l). (26) 



r 2 2wr 2 + 21 + 3 (2wr 2 + 21 + 3) 2 

The only difference between the equations ( 12~T]) and ( 1261) is that in the latter the 
parameter / is modified as / + 1 besides the constant term. As a consequence, from 
the known solutions, ()3]) and (@J, we can generate the solutions of the system fl25l) 
straightforwardly, that is 

-0+1), ' ' 



$ { nkr) = N n>l 



2^ + 2l + 3 L ^' [UJrle 



with 



and iV, 



£< 2 ) = £ n , m - w(2/ + 1) = Aco(n + 1) 

1/2 

= N n i + x is the normalization constant. 



ui l+ 2n\ 



{n + l + f)T(n + l + 2 



Now we investigate the action of SUSY operators, A^, on the eigenfunctions & n \ 



(27) 
(28) 



with the help of 
is 



4+cf> (1) 



For this purpose let us calculate the action of A + on $ 



nl+V tnat 



(2wr 2 + 2/ + 3) 
Evaluating the right hand side of ( 129|) we find 



■e 2 



(29) 



4+d> (1) 

A V n,/+1 



X 



alj n+l 



dr 

N n r +1 r l+2 



+ 2ur 



21 + 3 



4wr 



(2cjr 2 + 2/ + 3) 



(<+§) 

n+1 



-e 2" 



(2wr 2 + 2/ + 3)^ (30) 
By using the properties and recursion relations involving Xi-Laguerre polynomials the 

terms Mde the square bracket ean be replaced by -? ^ + 2 ' + ^< " ' = " + *' 



r V 2wr 2 + 2/ + 1 



sec 



Appendix A for details). Substituting this result in ( l30i) and simplifying the 



resultant expression we arrive at 



(31) 



Now let us evaluate the action of other operator A on the eigenfunctions that 



is 
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Here also while expanding the right hand side of ( I32p we find 



dr (2wr 2 + 2Z + 3) n+1 



1 , ,-2 



The terms inside the square bracket can be replaced by — 2ur ( 2^+21+3 ) (see 



(2ur 2 + 21 + 1) 



e—^ T . (33) 



Appendix A again for the details). With this identification, equation ( 1321) can be 
simplified to 

= - Ci^iW. (34) 

Equations fl3Tj) and ( l34"j) are the intertwining relations that relate the system ffT9l) 
and f |24l) . These intertwining operators while commuting yield = 2W'(r). We 

recall here that in the case of harmonic oscillator one finds that W'(r) = 1 and so 
the ladder operators associated with the Hamiltonian straightforwardly provide the 
Heisenberg-Weyl algebra. As a consequence the ladder operators of the harmonic 
oscillator perfectly act as annihilation and creation operators. However, in general, 
the factorization operators are not the ladder operators of the system for other than 
harmonic oscillator. The latter situation leads to the (nonlinear) polynomial algebras 
[36] . In particular, the super symmetric partners of the harmonic oscillator admit 
distorted versions of the Heisenberg algebra [37]. These nonlinear algebras can be 
linearized through the methods given in [38]. However, to the authors knowledge goes, 
for the shape invariant potentials other than harmonic oscillator, Balantekin algebraic 
method J26j|39] is more versatile to derive the necessary annihilation and creation 
operators. These operators frequently lead to either one of the following algebras: 
Heisenberg-Weyl algebra, SU(1, 1), SO(2) and (/-deformed algebra depending on the 
energy spectrum of the underlying shape invariant potential. As far as the present 
problem is concerned, on the Fock space, the commutation between the operators gives 
[A _ ,A + ]$ n ^ = 4uQ n j, which tells us that the SUSY operators A ± change both the 
integers n and I (vide equations ( I3T1) and (1341)). To capture perfect raising and lowering 
operators of the potential V\{r) (which should increase and decrease only the number 
of particles (n) by one and not the angular momentum (I) values) we proceed through 
the following way. 

2.2. Shape invariance property 

We attempt to explore the ladder operators by adopting Balentenkin's method |26j. In 
this approach, one can derive the ladder operators by analyzing the shape invariance 
property exhibited by the potential under investigation. The potential ([21!) can be 
shown to be shape invariant under the condition 

V 2 (r;l) = V 1 (r;l+ 1) +4u. (35) 

From this relation, we can define two new operators, namely B + and B_ , which are 
of the form 







B + = A + T{1) = A + e^ } B- = T~\l)A- =e~^A~ (36) 
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respectively. The commutation relation between these two new operators yields 

[B-,B+] = B_B + - B + B_ = e-^A-A + e^ - i f i. (37) 
Recalling the operators identity [4"0] . 

e A Be- A = B + [A,B] + ±[A, [A, B}} + (38) 

we can evaluate the first term appearing on the right side in equation ( 1371) . The result 
shows 

a - - , a d 2 99 1(1 + 1) 8u 
e~aiA~A + eai= — — 7 + u V + 



dr 2 r 2 (2ur 2 + 21 + 1) 

iy + i) 

(2cor 2 + 2l + l) K ' K ' 

On the otherhand evaluating A + A~~ we find 
i . i rf 2 99 Z(/ + l) 8a; 16w(2Z + l) /T w s 

Subtracting (j4"0|) from (|39|) . we get 

[S_,5+] = 4w. (41) 

The relation (14ip confirms that the operators I? + and £?_ form Heisenberg-Weyl 
algebra, which may also be observed from the energy spectrum of the potential Vi, given 
in ( 123]) . which is linear (vide equation (1231)). We note here that the new operators B + 
and B_ act as the ladder operators of the extended radial oscillator potential as in the 
case of harmonic oscillator potential [26] . 

From the above we can also establish 

£_<t>« = -Vi^ = -y/4u(n+l) (42) 

Further, redefining these two operators, 5_ and B + , in such a way that 

' B_, al = --^=B + , (43) 



we can show that 

M&g = d,,, ajsgj = v^TT $<V (44) 

The above relations confirm that these operators perfectly annihilate and create the 
eigenstate <3>^j by absorbing and emitting one photon. 

Using these ladder operators, one can construct coherent, squeezed and other type 
of states and analyze classical/non-classical properties exhibited by the radial oscillator 
potential (JT|). In the following we carry out these studies. 



3. Squeezed coherent states 



In this section, using the above ladder operators, we construct squeezed coherent 
states of the three dimensional isotonic oscillator. The squeezed coherent states of 
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the three dimensional spherically symmetric oscillator can be obtained by taking a 
tensor product of squeezed coherent states of radial excitation with squeezed angular 
momentum coherent states 130 1, that is 



|^a) = |^,a r )®||,a>. (45) 

The squeezed coherent states which generalize both the coherent and squeezed 
states are defined to be [27] 

\z,a) = S(z) D(a)\0) = D(a ) S(z)\0), (46) 



where S(z) = exp \{z*a 2 — za) ) is the squeezing operator with z = Re 1 ^ and 

D(a) = exp {acv — a*a) in which the operators a and a) represent annihilation and 
creation operators respectively. These operators satisfy the relations given in (I44p . 
Here a = «o cosh R + e 1 ^ sinh R, where R and 4> are the squeezing parameters and «o 
and «q are the coherent parameters respectively 

3.1. Radial part 

To begin with we evaluate the squeezed coherent states of the radial part, that is |£ r , a r ), 
by using the definition 

\£ r ,a r ) = D{a 0>r )S{z r )\0). (47) 

The squeezing operator can be disentangled as [27] 



1 2 
G(v ) tank Ry e^al — lncoshfi r aja r | tanhi? r e~*4ra% (48) 

cosh R r 



From ( 148)) . we find 



S(z r )\0) = I— tanhH r af| ^ (4g) 

cosh R r 

Substituting (T49|) in ( 147|) . we obtain 

|£r,Or) = -4^-^(oo,r) e^O), (50) 

cosh R r 

where we have defined £ r = — e^ r tanhi? r . Using the identity D^(a ^ r )D(a Qtr ) = 1, we 
can rewrite (!50|) of the form 

|£,,« r ) = D(a , r )e^ 2 f)^ a )D(a , r )\0). (51) 

cosh R r 

The combined action of D(a ,r) e^~ ar Z)'(a!o,r), yields 

J D(«o,r)e^ a ' 2 J D t (a 0ir ) = e^ {lll -*l, r ? ~ e 4pt 2 -< r?r at _ ^ 
Substituting this result in ( l5Tj) and simplifying the latter we obtain 



where a r = cosh i? r a , r + e^ r sinh i? r ckq r . 



Kr,«r) = I / = I /' 53 

cosh R r cosh it r 
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We can replace the exponential part in (153]) by Hermite functions, H n , n 
0,1,2,3,..., that is M\ 



e 



2° r + cosh ii r ar I V / / Sr I 



£ n! -f ) * < 54 » 

n=0 v 



so that equation (1531) now becomes 



oo H n 
n=0 ' ^ 



Qry/l-|€r| 2 



<Xry/l-\tir\ 2 



v 71 

?i=0 



2 # 

Thus the normalized squeezed coherent states of the radial part are found to be 



00 ^ n I ^/ _2g J / £ \ n / 2 

1^) = ^ V r / -| | n >, (56) 



Qr-y/l-lgrl 



n=0 



2 t 

where N r is the normalization constant whose exact value is given by 



N r = (l- |^| 2 ) 1/4 



exp 



~{a 2 r C + a?£ r + 2|a r | 2 ' 



, < |f r | < 1. (57) 



In the following, we construct squeezed coherent states of the angular part. 
3.2. Angular momentum part 

The angular momentum part of the Schrodinger equation (Q with d = 3 admits spherical 
harmonics Yi m (9,<p) as the solution which is defined to be simultaneous eigenstates of 
the operators L 2 and L z with the eigenvalues 1(1 + 1) and m, that is 

L 2 \l,m) = l(l + l)\l,m), (58) 
L z \l, m) = m\l, m), (59) 

where (8,4>\l,m) = Y^ m (9,<f) an d \l,m) are nothing but the angular momentum states. 

The angular momentum coherent states can be constructed by expressing the 
angular momentum coordinates, L + ,L_, L z and L 2 , in terms of the new annihilation 
and creation operators, that is a± and a± |!2"lPf3"] : 

(60) 

-(h + -h_), (61) 

(62) 





= a±a^, 




L z 


= g ( S + a - 


r — dLa_ 


L 2 


1 ,~ 
= ~ 2 (n + + 
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where h + + n_ = h is the number operator. The action of the operators a + and a_ on 
the number states is given by [42J 

a + \n) = y/n+ \n + — 1, a_ |n) = y / nr|ra+, n_ — 1), (63) 

a* + \n) = \/n + + 1 |n+ + 1, at |n) = a/ 77 - + 1 \n + , n_ + 1). (64) 

The operator L_ act on the states \l,m) as follows: 

L_\l,m) = y / (l + m)(l-m + l) \l,m-l). (65) 

The squeezed coherent states associated with the angular momentum part can be 
obtained by using the definition [29] . 

\i,a) = D(a o )S(i)\0,0). (66) 

Since we are dealing two modes, we consider the above equation is of the form 

|f, a) = D + (a 0j+ )^ + (e+)£ ) -K-)^-(e-)iO,0}, (67) 

with £± = — tanh R± e l<?i± and a± = «o,± cosh R± + «o ± sinh R± e 1 ^ are the squeezing 
parameters and displaced coherent parameters respectively. 

By adopting the calculations given in the radial part we find the combined action 
of displacement and squeezed operators produce (vide equations fl63l and ([54]) ) 



£>+(a ,+ )£+(£+) = ex P 



( a+ 



Vcoshi?, + 2 + 



(6* 



Recalling the identity, H n (x)— = e 2tx * 2 , we can express the operator appearing in 

n=0 

(I68p as an infinite series in Hermite polynomials, that is 



By induction, we find 



n+=0 



n + ! 



n+/2 



(69) 



£_(ao,_)&-(£-) = £ 



n_=0 



n_/2 



(70) 



Substituting (1691) and (1701) in (1671) . we obtain 

\la) = N ± ' 

n+,n_=0 



x 



n+/2 



"6 



n_/2 



a; aL |o,o), 



(71) 



where N± is the normalization constant which can be fixed as 

1 



iV± = (l-|e + | 2 ) 1/4 (l-|e-| 2 ) 1/4 exp 



-(a 2 + e + + alt + + 2\a + \ 2 ) 



x exp 



■-{a z X- + oi*_ f_ + 2|a. 
4 



0<|£±|<1, 



(72) 
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through the usual procedure. 
Using the identity, 

\l,m) 



it \l+m{rf \l- 



=|o,o), 



(73) 



y/{l + m)\{l-m)\ 

with the restriction n + = I + m and n_ = I — m, we can rewrite the squeezed angular 
momentum coherent states (1711) in the form 



\i,a) = N ± j2 E — 

z=o m=-Z 



x 



£+771 



-e- x 2 



V / (/ + m)!(/-m)! 

Z — m 

■e 



|Z,m). 



(74) 



Finally, substituting the squeezed coherent states of the radial part (vide equation 
( 156]) ) and the angular part (vide equation ( 174"]) ) in (145]) we obtain the squeezed coherent 
states of the three dimensional system ( 1T2]) which in turn reads 



oo oo I 



^U+J^S 



v 3 ^ 



l—m 



<*- \A-I€-I ; 



n=0 /=0 m=—l 

7 



^/n!(Z + m)!(Z — m ! 



|n, Z, m), 



o<ie r |<i, o<iai<i 

The normalization constant, A^ )Q ,, is given by 

4 

I2M1/4 



(75) 



[(i-ie,i 2 ) (i-ie+l 2 ) 

1 



x exp 



cxp 



exp 



— -(a££* + a* £ r + 2|a r 



--(alC + a_ £- + 2|a!- 



o<|£ r .|<l ! o<|ai<i- 

Equation (I75p can also be written in a more compact form, namely 



oo oo I 



|£, a) = N^ ta 22 zJ zJ CnCi,m\n,l,m), 



(76) 



(77) 



n=0 i=0 m=—l 



with 



(78a) 



and 



Q,r. 



v /(/ + m)!(/-m)! 



i + m 



(78b) 
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The squeezed coherent states given in (1771) are expressed in terms of both number 
states (\n)) and angular momentum states (\l,m)), that is the states are expressed in 
terms of bound state solutions i(j(r,9,(p). In other words a complete description of the 
three mode squeezed coherent states is established now in terms of three independent 
squeezing parameters (£ r , £+,£_) and three independent displaced coherent parameters 
(oi r , aq_, <x_). By varying these parameters we can analyze both the classical and non- 
classical nature of these states. We recall here that whenever the quantum states exhibit 
Poissonian statistics they said to possess classical nature |44j. Any deviation from 
this behaviour represents the non-classical nature of the states. To proceed further we 
restrict the range of values of £ 3 - to be < |^| < 1 and the range ao,/s be — oo to oo 
since normalization constant N^ a is defined in the range < < 1, where j = r, +, — . 

In the following section, we analyze the classical/non-classical nature of the 
obtained squeezed coherent states. 



3.3. Completeness condition 



In this sub-section, we demonstrate that the three mode squeezed coherent states 
([77]) satisfy the completeness condition. The three mode squeezed coherent states are 
represented in terms of three coherent parameters (oc jS, j = r, +, — ) and three squeezed 
parameters (£js, j = r, +, — ). Now we prove that the squeezed coherent states resolve 
the identity operator [15] 

(79) 



J «| d a = I, 

when the integration is carried over the entire space of a (a r , a+, oj_). In this analysis 
we consider a varies only with respect to a 0j? s and treat £-s are arbitrary constants. 
For the squeezed coherent parameters, 



O ; 



we can show that 



dajda* 



a o,j - 




0- 


- IOI 2 


dasj 








da* 


8a* 







, j = r, 0, 0, 



daojdaQj = daajda^j. 



10) 



51) 



Substituting the expression given in (1771) in the integral ([791 . we get 

/oo oo !' oo oo I 
\£,a)(£,a\d 2 a=^2Y^ Yl Yl Yl Yl \ n ' ' l> > m ')( n ' ^ m \ 

n'=0 l'=0 m'=-V n=0 1=0 m=-l 




c* n ,c* lrnl c n ci, m Nl a d 2 a r d 2 a + d 2 a_. (I 
To start with we evaluate the radial part (which we call as Gr) in ([82]) . that is 



■2) 



oo oo „ 

Gr = J2Y\ n '^ U \ / C *n> C n N ?J 



71=0 n'=0 
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tt ^a^(-f f (4)7- h' 2 - 1 

n=0 n'=0 



X # n , 



-2£ ; ~" v 

Substituting ( 180|) and (IHTj) in (!83|) and separating the integral (which we call as ii) from 
the summation part, we get 

, 1 




exp 



-\ a o,; 



£ r O! 0>r . + Cr a 0,r 



^O.r / C , tt 

+ \/-y«o,r I ii n 



+ \/-7r a or I da 0ir da* 0r . (84) 



To evaluate the integral (J84J) we introduce the transformation 



J ;/J " Jl ^ and 



<j v) + My 3 



so that 



i-IOI 2 

dapj dap j 



-26 1-I6f 



daojda j 



dyj dy* 
9a h da o,j 



dyj dy* 

In the new variables the integral f )84|) reads 



-2|£-| 



55) 



(86) 



-2ie r 




exp 



ie,i 2 (^ 2 + i/; 2 ) + 2^ik 
i-ie,i 2 



H n {y r )H n ,{y* r )dy r dy;. (87) 



With another change of variable, z = z -, ^ y r} the integral ( IH71) can be brought 

\/l-|fr| 2 

to the form 



-2z 



1 " % 



exp 



x#„ 



^n[2/*] / eX P 



V 1 - ICr 



16-1 



-2T + 2z 



dz dy*. 



With the identity [H 



a 



(89) 



the second integral in flHHl) can be evaluated and given in terms of Hermite polynomials, 
that is, 



h 



-% 



exp 



1-krl 2 



fl" n [-y r *]^[y r *]dy r *. (90) 



Now using the orthogonality property of the Hermite polynomials, 

J e' x2 H n [x]H' n [x]dx = y/n 2 n n\ 5 n>n ,, 



(91) 
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where 5 n ^ n i is Kronecker-delta function, we find 

n 



h = , 2 " = (~ ) n\5 nn ,. (92) 



Substituting this expression, (192]) . in Gr (vide equation (|83|) ) and simplifying the 
resultant expression, we get 



n=0 n'=0 N 7 

From the above we observe that radial part of the integral provides 

oo 

G R = -2m^\n)(n\. (94) 

n=0 

In a similar way we can evaluate the integrals involving the other two variables, 
namely 9 and (since the procedure is exactly the same as that of radial part we do 
not repeat the details here). Our result shows that 

OO V OO / /././. 

G A = ^2 ^2 ^2^2 \l'> m ')(h m \ / / / c* n ,c* llm ,c n ci )m Nld 2 a r d 2 a + d 2 a^ 

l'=0 m'=-V 1=0 m=-l J J J 

oo I 

= -4:7v 2 Y^ \ l , m ){ l M- ( 95 ) 

1=0 m=-l 

The integral ( 18 2 p is completely evaluated now. The resultant value turns out to be 
G = Gr x Ga = 'S^ 3 Er=o Ei^o EL=-i \l,m)\n)(n\(l,m\. This in turn confirms that 
the states |£, a) resolve the identity operator, that is 

oo oo I 

\n,l,m)(n,l,m\ = I. (96) 

n=0 1=0 m=—l 

The result ensures that the constructed three mode squeezed coherent states form a 
complete set. 



4. Non-classical properties 
4-1. Mandel's Q parameter 

In this sub-section, we study Mandel's Q parameter for the three mode squeezed coherent 
states (I77|) . We evaluate Mandel's Q parameter associated with each mode |4B] 

Qj = tt4 - (fij) - 1, j = r,+,-. (97) 
\ n j) 

Here, rij denotes the number of particles in the respective mode which can be uniquely 
determined from their associated number operators hj = cijCij. The action of the ladder 
operators a] and cij on the number states \rij) are given in equations ( 1441) . ( )63l) and (1641) 
respectively. 
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To evaluate Mandel's Q parameter we consider the squeezed coherent states be of 
the form 

oo oo oo 

|£, a) = N^ a ^2 ^2 E c n c n+,njn,n + ,n^), (98) 

n=0 n+=0 n_=0 

where the constants c n and c n+i „_ are given in (178a[) and ( 178bl) with Z ± m should be 
replaced by n±. To find the expectation values (fij) and (n 2 ), as we did earlier, we first 
calculate the radial part. Doing so we find 



oo oo oo 



{n 



* * 

C n C n + ,n- C n- c n + ,n- n 



and 



•)=■«?,. EE E 

n=0 n_|_=0 n_=0 

= (|^ [2 1 ,!) [l^l^ 1 + l^l 2 ) + ZrCtf + + \£r\ 2 

oo oo oo 



c ra c n+,n_ c n c n + ,n_ re 



n=0 n+ =0 n_ =0 



*+2 



(1 _ ^ |2)2 [Kl 4 (i + 1^| 2 ) 2 + |£ r | 2 (2 + |£ r | 2 ) + (£> 2 + HrK , 

+(2(1 + |a r | 2 )(l + \i T \ 2 ) + 2\i r \ 2 W r a 2 r + £ r a? ) + |a r | 2 (l + 8|£ r | 2 + 3|£ r | 



(99) 
(100) 

(101) 
(102) 



The details of evaluating both the expressions fl99|) and (110 II) are given in the 



Appendix B 



Since n±|n + ,n_) = n±|n+,n_), the expectation values (n±) and (n±) associated 
with the angular part also provide the same expressions given in f II 1) and (11021) . Hence, 
in general, we can write the expectation values (hj) and (n 2 ) are of the form 



(flj) 

(n 2 ) 



^I-IOI 2 ) 



(1-I0I 2 ) 2 



\<*j\ 2 (l + \£j\ 2 ) + Zj<xf +$a? + |&| 2 



(103) 



«i! 4 (l + l^ff + fcf(2 + l^l 2 ) + (e*« 2 + ) 2 



+(2(1 + | aj f)(i + ion + %r)(ex ) + M + %i + %i 4 )j ,(104) 

«o,j-0«o,i , . 
where a,- = — . and 7 = r, +, — . 

3 V^W 2 

From the expressions f 1 1 U 3 j) and (11041) we can analyze the Poissonian (Qj = 0), sub- 
Poissonian (Qj < 0) and super-Poissonian (Qj > 0) nature of the states in each mode. 
We can also calculate the parameter Qj of the three mode squeezed coherent states which 
basically depends on the parameters aoj and £j. The obtained numerical results are 
depicted in figured! The figure shows that when \£j\ < 1 and aoj = 3, the (3-P arame t er 
takes positive and negative values which in turn confirm the super-Poissonian (Qj > 0) 
) and sub-Poissonian (Qj < 0) nature of the states. The sub-Poissonian statistics 
indicates the non-classical nature exhibited by the states. 
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-0.5 



0.0 0.2 0.4 0.6 0.3 

If j I 

Figure 1. The plot of Mandel's Q parameter for each mode of the states with Oq j = 3. 



4-2. Quadrature squeezing 

The squeezed coherent states given in (177)) are expressed in terms of the variables r, 9 
and (p. As far as the radial part is concerned one can analyze the squeezing in the 
observables such as generalized position (w r ) and its conjugate momentum (p r ). To 
show the squeezing in the angular variables 9 and we have to consider the angular 
momentum quantities L x ,L y and L z . In the following, we analyze the squeezing in 
position and momentum coordinates which are defined as [47] 



-j=(al + a r ), 



Pr 



V2 



(105) 



To analyze the squeezing in the quadratures w r and p r in which the Heisenberg 
uncertainty relation holds, (Aw r ) 2 (Ap r ) 2 > |, where Aif5 r and Ap r denote uncertainties 
in w r and p r respectively, we introduce the following inequalities, that is 

h = (a?) + (af ) - (a r ) 2 - (a\) 2 - 2(a r ){a\) + 2{a\.a r ) < 0, (106) 
h = ~ (al) - (af) + (a r ) 2 + (al) 2 - 2<a P )(aJ) + 2(ala r ) < 0, (107) 

which can be derived from the squeezing condition (Af) 2 < | or (Ap r ) 2 < | 
by implementing the expressions given in (110511 . The expectation values should be 
calculated with respect to the squeezed coherent states |£, a) in which the squeezing 
property has to be examined. 

For the squeezed coherent states (1771) . we obtain the following values for the 
quantities which appear in the equations ( 11061) and (1107j) . that is 



(a r ) 



( OL r + |Cr| 




, (al) = 


I a* + \£ r \ 




w 1 - 16 






IV 



(ala r ) 
where a r = 



(af) 



(i-ie,i 2 ) 



\Zr\ 2 ' XW l-|£r| 2 ' 

(|a r | 2 (l + %\ 2 ) + %\ 2 + a% + afc) 



(108) 

(109) 
(110) 
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With the expressions given in fll08j) - fll 10p we evaluate the inequalities ( 1106)) and 
( 1107P numerically and plot the outcome in figure [2] with a ,r = \<^o,r\€ l9 °' r ■ From figure 
|2] we observe that the identities given in equations (1106!) an d (11 07ft for the squeezed 
coherent states |£, a), satisfying the uncertainty relation, shows I 2 < and l\ > 0. This 
in turn confirms the squeezing in the quadrature p r , for all values of atQ >r . 

As far as the angular part is concerned the squeezing in the quadratures L x or L y 
can be analyzed through either one of the normalized quantities Sl x < or Sl v < 

where SV = 2A ^-~^ Z ^ and Sr = 2AL ,% rr^ 1^81 - We can express the operators L x 

and in terms of L + and L_, namely 

4 = i(£+ + L_), 4 = (in) 

so that the uncertainty of L x and L y can now be expressed in terms of L + and L_, that 
is 

= j {(L\) + (Li) + (L + L-> + (L-L+) + + <L_» 2 ) , (112) 

(A^) 2 = (L 2 y ) - (4) 2 

= - \ ((Ll) + ) - (L+L.) - (LM) ~ ((L+) - (L_)) 2 ) . (113) 

For the squeezed coherent states (1771) . the expectation values in (I112p and (1113)) are 
found to be 
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5b s 





3 4>0.S 



(b) 



Figure 3. The plot of normalized quantities (a) Sl^ and (b)Sx with = 0.1 

and oto = 1.3. 



L 



L. 



+ 1 , (116) 




((<*+ - C+a* + ) («; - e+a+) + 



x ((a_ - (al - £a_) + IC-I) 



2 



(120) 



where a± . 

We evaluate the normalized quantities Sl x and Sl v numerically by substituting 
the expectation values (1114j) - (1120l) in the uncertainties (11121) and f l 1 1 31) . We plot 
the numerical results in figure |3] where we have considered ao,+ = x + + iy + with 
£ + = = 0.1 and ao,_ = 1.3. For this choice of parameters figure [3] shows that 
Sl x > and Sl v < which explicitly demonstrates the squeezing in L y . 



5. Wigner function for 3-dimensional isotonic oscillator 

In this sub-section, we evaluate Wigner function (W(x,p)) of the squeezed coherent 
states ( 1751) . The Wigner function (a quasi-probability distribution function), which 
was introduced as quantum corrections in classical statistical mechanics, normally takes 
negative values in certain domains of phase space so that it cannot be interpreted as a 
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classical distribution function which is non-negative by necessity [15] . 

The Wigner function for the single mode state is described by the density operator 
p which can be written as |31j 

W(() = 2 Tr[pf (C, 0)], T(C, 0) = D(C CV^^CC, O, (121) 

where D((, (*) and p are the displacement and density operators respectively and T is 
the complex Fourier transform of the s-parameterized displacement operator, D(i], s), in 
which rj is the coherent eigenvalue with s = corresponding to the Wigner distribution 
function. 

Since p = |£, a)(£, a|, the Wigner function for a single mode can be explicitly 
calculated from 

W(() = 2(Z,a\T((,0)\a,0. (122) 
Prolonging the definition (11 22ft to three orthogonal modes [32], we get 

3 

W({d, C 2 , Cs}) = 8 Tr[p HfiiCi, 0)], fi(ti, 0) = D{Ci)e l ^D^iQ), (123) 

i=i 

where Ci = Cn C2 = C+ an d C3 = C— By applying this definition, (11231) . to the three 
mode squeezed coherent states fl75|) . we obtain 

oo oo Z V 

wac,( + ,c-}) = Ni a e E E E 4'C„4 im ,Q, m 

n,n'=0 Z,Z'=0 m=-Z m'=-2' 

x (Z', m'| (n'|f (C,, 0)f (C+, 0)f (C_, 0)|n)|Z, m> , (124) 

where c n and q >m are given in fl78aj) and (178bl) . The above equation can be evaluated 
by separating the radial part from angular part, that is 

00 00 1 v 

W({( r ,( + ,(-}) = Nl a E E E E <,c n cl m ,c l>m (n>\f(( r ,0)\n) 

n,n'=0 l,l'=0 m=-l m=-V 

x(l',m'\f(C + ,0)f(C-,0)\l,m). (125) 

The operation of the dual state (n'| on the state T(£ r ,0)\n) may be known as 
transition probability, expressed in terms of n, has already been reported in Ref. [3 1 J . 
The result shows that 

(n'|T(C,,0)|n) = e- 2 K"l 2 ffl^ 2"-"' +1 (-l)"'(C) n - ri '^r"'(4|C| 2 ), (126) 

where L™~ n ' is the associated Laguerre polynomial of degree n. In a similar way, the 
quantities corresponding to the operators T(£ + ,0) and T(£_,0) can also be found in 
terms of the number states n±. The resultant expressions are turned out to be 

(n' + |f(C + ,0)K> = e" 2 ^l 2 ^y /2 2« + -«;+ 1 (-l)«;( C ;)"+-^L:r n ' + (4|C r r), (127) 
and 

<nl|f(C-,0)|n_) = e- 2 l«-l 2 ^y / % w --"- +1 (-l)"'-(C) n --^C~ n '-(4|C-r). (128) 
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Substituting (1 126ft - ( lT2"gj) in f )125p . we can obtain Wigner function for the squeezed 
coherent states f!75|) with n + = I + m, n' + = I' + m! and n_ = I — m, n'_ = V — m! 
in the form 

oo oo I V 

^({c,c + ,c-}) = e-i^i 2 -i^i 2 -i^i 2 iv|, Q E E E E 4, Cn 4 im , Q>m 2^'- 2 ^ 

n,n'=0 M'=0 m=-l m'=-V 

' n'\(V + m')\(V - m')\\ l/2 



n\(l + m)\(l — m)\ 

|2\ Tl+m-V-m 



n' I /-*\n—n' / /■* \ l+m—l'—rn' 



(-if (C)"- n (CI) 



(C 



xL»r B ^4|Cr| a )LK l '- m X4|c + r)4-V 4W (4|C 



(129) 



where = ^ + Pj, & = — f tanh -Rj-e^ and <x,- = a j cosh Ri + ag^e*^' sinh Rj, j = r, ±. 
On evaluating the right hand side (the details are given in the Appendix C), we 



find 

^({C„C + ,C-}) = nLex P 



+ 



1-16 1 2 



- 2 k 



(130) 



We analyze the Wigner function graphically by considering £ r = C+ — C-> 
£ r = £ + = and «o,r = «o,+ = «o,-- With this assumption, the Wigner function 
of the three mode states can be seen as one mode state. We then calculate the function 
W(x,p) numerically and plot the outcome in figure H] for a specific value of R, (ft and «o 
with C = x + i p. 




Figure 4. The plot of Wigner function corresponding to squeezed coherent states (|75|) 
with (a) R — 0.7, = and ao = 0.5 

The plot confirms that the Wigner function takes smooth shape in which the 
squeezing effect can be seen. 

6. Conclusion 

In this paper, we have constructed squeezed coherent states and studied their 
nonclassical properties associated with the 3-dimensional generalized isotonic oscillator. 
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Since the potential under consideration is a spherically symmetric one, we dealt 
the radial and angular momentum parts separately As far as the radial part is 
concerned we observed that it is nothing but the newly found extended radial oscillator 
whose eigenfunctions are expressed in terms of the recently discovered Xx-Laguerre 
polynomials. We have shown that this radial part exhibits shape invariance property. 
Using this property and implementing Balantekin's method, we have obtained two new 
operators which in turn constitute the intertwining supersymmetric operators. We 
have demonstrated that these operators produce Heisenberg-Weyl algebra and perfectly 
annihilate and create the eigenstates of the radial part. Using these ladder operators 
we have constructed squeezed coherent states of the radial part. As far as the angular 
part is concerned we have used Schwinger's representation to define the creation and 
annihilation operators. We have derived the associated squeezed coherent states by 
considering the disentangled form of two mode squeezing and displacement operators 
with certain restrictions on the allowed values of angular momentum variables I and m. 
Finally, we have expressed the three mode squeezed coherent states of the generalized 
isotonic oscillator as a tensor product of the squeezed coherent states of the radial part 
and squeezed coherent states of the angular momentum part. To make our results 
more rigorous we have proved that these three mode squeezed coherent states resolve 
the identity operator. We have explicitly illustrated the squeezing properties of the 
constructed 3-dimensional squeezed states both in the radial part and in the angular 
part. As far as the radial part is concerned we have considered the conjugate variables, 
generalized position w r and the momentum p r and proved the quadrature squeezing. 
Similarly, to show the squeezing in the angular momentum we have considered two 
normalized quantities, namely SV = 2 ( AL ^ zlihli anc i 57 = 2 ( AL ») ~K Lz )l anc j 
demonstrated explicitly the squeezing properties possessed by these states in the angular 
variables. In addition to the above, we have evaluated Wigner function of these squeezed 
coherent states and shown that these states exhibit squeezing property. Since the ladder 
operators of this 3-dimensional system act on the eigenstates only linearly, one can 
deform these operators nonlinearly and construct certain non-classical states such as 
nonlinear coherent states, nonlinear squeezed states, generalized intelligent states and 
so on. The details will be presented elsewhere. 
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Appendix A. Evaluation of ((30j) and ( 155ft 

To evaluate the terms inside the square bracket in (|30|) we recall the following property 
associated with the X\- Laguerre polynomials [5]. The differential equation ([5]) can be 
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factorized by two operators [5], namely 

(x + k + 1) 2 d 



and 



Mv) 



B k {y) 



(x + k) dx \x + k + 1 



-(2/' - y) + ky. 



(A.l) 



(A.2) 



(x + fc + l) 

One can unambiguously check A k B k (y) factorizes the equation flSJ as A k B k (y) = X(y). 
Rescaling the variable x = uir 2 and the constant k = I + \ in (LA. 21) . we get 



Bi(y) 



2ur 2 + 21 + 1 
2 \ 2ur 2 + 21 + 3 



d i + | 2,r 

dr 



21 + 3 



Using the identity [5] BiL 
brought to the form 



(2ur 2 + 21 + 1) 



(A.3) 



nL„ +1 with y = , equation ( ]A.3j) can be 



aLl n+l 

dr 



+ 2ur - 



21 + 3 



+ 



Auir 



(2ur 2 + 21 + 3) 



L 



n+l 



2 /2wr 2 + 2Z + 3 
~ ~r \2ur 2 + 2l + 1 

x(n+l)L^p. (A.4) 



We use this expression to simplify the equation ( 130]) . 

To evaluate (13"3"|) . we again rescale the relation OA.ip in such a way that x = ur 2 
and = I + \ so that the identity (jA.lj) now becomes 



Mv) 



2ur 2 + 21 + 3 



2ur V 2ur 2 + 21 + 1 
i+h) 



Aur 



dy_ 

dr (2ur 2 + 21 + 3) ' 



(A.5) 



Recalling another identity A[L n+1 
rewritten as 



~(l+2.) 

Ln 2 with y = L\ +l , equation (]A.5j) can be 



n+l 



4wr 



(/?' 



(2wr 2 + 2/ + 3) 



>+*: 

n+l 



—2ur 



f 2ur 2 + 2l + l \ -(/+§) 
\2o;r 2 + 2/ + 3 



-^n 



(A.6) 



This relation is used to simplify the intertwining relation ( 133]) . 



Appendix B. Evaluation of ((991) and (TIuTTl 

In the following, we discuss the method of evaluating the series (n r ) given in ( 1991) . that 
is 

oo oo oo 
n=0 n+=0 n_=0 

With the definition A^ 2 = X!„%o ^n°_=o c n+,n_ c n+,n_, where A| is given in ([72]), the 
above equation (199]) can be written in the following compact form 

oo 

(h r ) = N 2 Y,c* n c n n. (B.l) 

n=0 
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Substituting (178aj) in (IB.li and redefining the summation appropriately, we obtain 



n, 



^ _ jy2 ^ H s+1 (x* r )H s+1 (x r ) / -|g r [ \ 

s=o '' V 2 / 



a+1 



(B.2) 



Recalling the integral representation associated with the Hermite polynomials, 
namely H n (z) = ^ J^ 00 e~* 2 (2; + it) n dt |JT], the above equation f ]B.2j) can be brought 
to the form 



(n r ) 



7T 



OO /"OO 



e Zl " ' 



OO ./ — OO 

OO 



X 



s=0 



2 (x r + izi)(x* + zz 2 ) 
-2|^ r |(ar r + wi)(a;*+i3!2)) 



SI 



dz\dz2- 



(B-3) 



We observe that the summation inside the integral can also be written as an 
exponential function. The non-exponential terms left over inside the integrals, that 
is (x r + izi)(x* + iz 2 ) can also be brought into exponential form through the identity 

(x r + iz{)(x* + %z%) = { ^ [exp \p(x r + izi){x* + ^2)]]} • As a result equation (1B.3|) 

now becomes 



p=0 



7T 



x - / / exp - z| + (p - 2|^ r |)(av + wi)(sj + iz 2 )] dzid.Z2 f (B.4) 
I dp \_J — 00 7-oo J J p= o 

To evaluate (1B.4[) we introduce the transformation 77 = p — 2|£ r |. In the new variable 
equation ( 1B.4|) reads 

-2|£ P |JV? f d 



(n r ) 



OO POO 



exp [— «x ~~ z 2 + ^(^r + ^l)(^* + iz 2 )]dzidZ2 



OO ./ — OO 



7r ^ dry 

To begin with we evaluate the following integral in ( 1B.5|) . namely 

/OO 
exp \—z\ + rj(x r + z^i)(x* + iz 2 )]dzt. 
-00 



(B.5) 



7?=-2|? r | 



(B.6) 



We noted that this integral is nothing but the Fourier transform of an exponential 
function. With this identification we find 

.2 



h 



7T 



exp 



2 * 2 2 

— £ 2 ^ «Y X ''' Z2 + ^l Xr l +«^r^2 



Substituting ( IB. 71) in the second integral in ( 1B.5[) . we get 

/>oo 

h= I exp [-z 2 y \l 1 dz2- 



(B.7) 



(B.8) 



The above integral also turns out to be the Fourier transform of an exponential function. 
The explicit integration leads us to 
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Substituting (IB.9P in (1B.5P and simplifying the resultant expression we arrive at 



(h r ) = -2\£ r \N? { ± 



: exp 



1 - 



ry-i I rp I 2 / ™2 _|_ ™ * 

'/I m 4 v^r' ~ r 

1 _ nl 



. (B.10) 



r,=-2|£ r | 



Now carrying out the differentiation and substituting the expressions 



and its complex conjugate in the resultant equation we obtain the 



following expression for (n r ), namely 

1 



(n r ) 



«r| (1 + 1^1) + « C + °C + \£ 



:i-i^i 2 ) 

We use this expression to evaluate the Mandel's Q parameter. 

Now we evaluate the expectation value {h 2 ) given in (llOip . that is 



(HOOD 



oo oo oo 



(rr 



n 2 . 



(B.ll) 



n=0 n+=0 n_ =0 



n_=0 u n + ,n_ <-ra + ,n_ , 



Separating the radial part from angular part by recalling N± 2 = =o S 
with N± defined in (172|) . the triple sum in (IB. lip can be reduced to the form 

oo 

(n 2 r ) = N 2 Y,<c n n\ (B.12) 

n=0 

Substituting f)78ap in flB.12|) and redefining the summation appropriately, we find 



(rr 



H s+2 (x*.)H s+2 (x r ) f-\C, 



o! \ 2 



3=0 



H j+1 (x* r )H j+1 (x r 



■\Cr 



3+1 



(B.13) 



The second term in (1B.13P is nothing but (h r ) (vide equation (1B.2[) ) which has 
already been evaluated (vide equation ( llOOp ). As a consequence, we confine our attention 
only on the first term in equation (IB. 13|) (which we call as G), namely 



J H s+2 (x* r )H s+2 (x r ) f-\£ r x ~ + ' 2 



Si 



(B.14) 



Using the integral representation, H n (z) = ^ e f2 (z + it) n dt , we can rewrite 
(1RT4D as 



G 



2 I C 1 2 roo pco 



71 



e Zl z *(x r + izxfixl + iz 2 f 



— oo J — oo 

OO / 



i-2|^|(^r + iz 1 )(x* + iz 2 )Y d ^ 



s=0 



Si 



(B.15) 
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By recognizing the sum appearing in ( IB.15P is nothing but an exponential function 
and the terms (x r + izi) 2 (x* + iz 2 ) 2 can also be rewritten as an exponential function we 
can bring equation to the form 



G 



TT 



( d 2 r f°° f°° 1 1 

x \l~2 / / exp[-z 2 - z% + (p -2\ti r \)(x r + izi)(x r + iz 2 )]dzidz 2 > .(B.16) 

To evaluate the integrals we introduce the transformation rj = p — 2|£ r | so that the 
equation (IB. 161) in the new variable reads 

4iV r 2 |£ r | 2 



G 



71 



( d 2 r f°° f°° ii 

I / exp i~ z i ~ 2 2 + V( x r + izi)(x* + iz 2 )]dz 1 dz 2 > . (B.17) 



Evaluating the integrals with the help of Fourier transform, we find 



G 



7T 



drj 2 



cxp 



7]\X r 



4 V r ' r 



1 



(B.18) 



V=-2\£ r \ 



Now differentiating the terms inside the square bracket two times with respect to 
r) and then substituting the expression for N r in the resultant equation and simplifying 
the latter we arrive at 

G = I m2 [k| 4 (l + m 2 ) 2 + (a 2 C + aftr) 2 + 4|a r | 2 |e,| 2 (2 + |^| 2 ) 

+ |£ r | 2 (l + 2|£ r | 2 ) + ((2|« r | 2 (l + |£ r | 2 ) + 1 + m 2 )(» 2 Cr + uftr)] ■ (B.19) 

Substituting the equations (IB. 191) and ( 11001) in (IB. 131) . we get the following 
expression for (n 2 ), that is 
1 

(l-l^l 2 ) 2 

+ (2(1 + |a r | 2 )(l + |£ r | 2 ) + m\ 2 )(Ca 2 r + ZrOL?) + \a r \ 2 (l + 8|£ r | 2 + 3|£ r | 4 ) .JM} 



(n 2 ) 



«r| 4 (l + m 2 ) 2 ) + |6| 2 (2 + |6| 2 ) + {C<4 + &•« 



To numerically evaluate the Mandel's Q parameter we use the expression ( llOip . 



Appendix C. Evaluation of Wigner function ( 11291) 



Let us first evaluate the radial part in ( 11290 . Separating the radial part from angular 
part we get 



W r {Q = e~ 2 ^ 2 N 2 E Hn(Xr) ^ AK) (2CT-' (- 

n=0 n'=0 



n/2 



c 

2 



n'/2 



I 2 ) 



:-ir 

(C.l) 



where x r 



^ry/l-l^l 2 
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To start with we evaluate the following sum appearing in (IC.ip . that is 

n/2 

LlT n '(4|C| 2 ). (C.2) 



OO TT I 



n=0 



III 



"6 



Recalling the integral and summation representation of the Hermite and the 
the associated Laguerre polynomials, namely H n (z) = ^= / 00 °°e~' 2 (2; + it) n dt and 

Yl^Lo L " = eW ^n( z ~ w ) w e can reduce equation (1C.2|) to the form 



7T 



x / exp [-t 2 + 2(x r + it)C ^f^^Xl K?' ( 4 Kr| 2 - 2y^2l r C(x r + it)) dt.(C.3) 

This expression can further be simplified by recalling yet another identity, that is 
L~ n (x) = ^j—. With this simplification, equation (10.3j> reshapes into 



Gi 



(-1)" 

r-OO 



exp -t 2 + 2(x r + it)Gy/-2Z r U\( r \ 2 -2 y^2YrC(x r + it)) dt. (C.4) 

-oo 



The second term in (10.41) can also be written as an exponential function through the 



relation (4|( r | 2 - 2^ = 2^(*(x r + it)) 
By doing so we find 

'-1Y' { d n ' 



£ [exp [p(4|C r | 2 - 2(x r + it)V=2frC)]] } 



p=0 



G 1 



n'\\pH \dp n ' 



exp 



4p|C r | 2 + 2(l-p)V = 2^C 



x / exp 



-t 2 + 2l(l-p)y^2j r C r t 



dt 



(C.5) 



p=0 



Now evaluating the integral using the relation e twx e x " dx = y/ne * , we obtain 



n'\ \ dp n ' 



exp 



4p|C,| 2 + 2(1 - p) v^Or + 2(1 - pfUl 



(C.6) 



p=0 



To proceed further we introduce the transformation, 8 = y/—2£ r Q(l —p), in ( IC.6I) . 
In this new variable 9, equation ( 1C.6I) reads 



n'\ 



exp 



4|Cr| 2 + X r + J—C 



-2 



x 



d n ' 

d(F 



exp 



9-X r - \l -^Cr 



(C.7) 



Expressing the differential part appearig in ( 1C.7I) in terms of Hermite polynomials 
through Rodrigues formula, H n (x) = (— l) n e x2 4-^e~ x2 [ 1 1] - we find 



G\ = n exp 



TV 



2U*r + 2 V /Z 2^,C 



(C.i 
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Substituting (10.8|) in the Wigner function (10. ip . we get 

H n ,{x*)H n ,{9) (-\i r 



W r = Ac exp 



n'=0 



TV ! 



(C.9) 



where 9 = >/=5£C - 



The summation in ( 1C9f) can be evaluated through the identity [2] 



£ 

n=0 



H n {x)H n {y)w r 



m 



y/l - 4u> 2 



exp 



Aw 2 (x 2 + y 2 ) — 4ira?/ 



4w 2 - 1 



H<2- ( ai °) 



On replacing this result in (10.9 j) . we obtain 
W r = N 2 exp 



-2ec 2 + 2v^;c-2iCri 2 + 

2 2 



-2 



-I (x 2 + x; 2 + 4|c 

2|Cr| _/ 
1 ^ 



I6| 2 - 



(on; 



The final task is to substitute the normalization constant, (137)1 . in ( 10. lip and 

ar a/1 — if r p 

simplify the resultant expression with x r = r ryur ~- On completing this task we 



arrive at the following expression for the radial part, that is 



W r = exp 



-2|Cr| 2 (|6 


2 + i) + 2(e r c 2 + cc 2 ) 


i - & 

2(C«rCr + ^<0" 


2 






2 



- 2 a r r + 



2 ^(^tCt ^rCr 



(C.12) 



In a similar way one can also evaluate the other two modes present in the Wigner 
function. Our result shows that (since the procedure is repetitive, in the following, we 
give only the final form of the expression) 



W± = exp 



-2ic±r(i£±i 2 +i)+2(£±e r+sc 



-2|a±| 



1-16 



(C.13) 



As a result, we obtain the following form for the Wigner function for the three 
dimensional generalized isotonic oscillator, that is 

W({( r ,( + ,(-}) = W r xW ± 



ITj =1 exp 



-2ioi 2 (ioi 2 + i) + 2fec; 2 + eKI) 



2N 



(C.14) 



We use this result to investigate the Wigner function numerically. 
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